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Two techniques of measurements of thin film magnetostriction are compared: direct, when changes 
of the substrate curvature caused by the film magnetization are controlled, and inverse ("indirect"), 
when the modification of the magnetic anisotropy induced by the substrate deformation (usually 
bending) is measured. We demonstrate how both the elastic strength of the substrate and the 
effective magneto-mechanical coupling between the substrate deformation and magnetic anisotropy 
of the film depend on different conditions of bending. Equations to be used for magnetostriction value 
determination in typical cases are given and critical parameters for the corresponding approximations 
are identified. 

PACS numbers: 46.70.De ; 62.20.dj ;75.70.Ak;75.80.+q ;85.85.+j 
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I. INTRODUCTION 

There are two alternatives for magnetostriction mea- 
surements of bulk ferromagnetic materials: it can be 
measured either directly, as it was initially observed from 
deformation under the magnetic field, or by an inverse 
way measuring the magnetic anisotropy induced by me- 
chanical stress applied to the sample 3 . In thin films 
the magnetostriction deformations are hindered by the 
substrate and only a secondary much smaller substrate 
deformation can be detected. By analogy to the bulk 
case, the direct technique uses the substrate deforma- 
tion caused by the modification of the film magnetiza- 
tion. The inverse technique (often called " indirect" , the 
term we find not very adequate) implies studies of the 
thin film anisotropy changes caused by the deformation 
of the substrate (see^i and references therein). 

Description of these techniques involves elastic prop- 
erties of both the thin film and the substrate and, in 
particular for the direct technique, it took a long time 
to find the right equation to extract the magnetostric- 



tion parameter; 
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. The entire mechanical prob- 



lem is so complex that even after the right solution 
was well established-- other contradicting theories are 
published 12 .. The theory for the inverse method ap- 
pears to be simpler because the shape of the substrate 
is controlled directly and one can consider only the thin 
film properties. Nevertheless, different expressions of the 
stress induced anisotropy can be found in different pa- 
pers without a clear explanation of their origin ^ 18 ! 21 . 
The difficulties come from an interplay between longitu- 
dinal and transverse deformations in response to uniaxial 
stresses related by Poisson's ratio. The problem is cru- 
cial for applications of cantilevers as sensors and actua- 
tors where the role of the Poisson's ratio was specially 
investigated^. 

The rigorous solutions for the bimorphous plate bend- 
ing includes so heavy and not transparent equations that 
for scientists of the magnetic community it is difficult to 



understand the reason why the Poisson's ratios of the 
substrate and of the film appear in so different forms. At 
the same time the direct experimental verification of the 
theory is not convincing due to the fact that the param- 
eters of thin films can vary significantly compared to the 
bulk onesi£. 

In this paper without repeating the complete deriva- 
tion we present simple criteria for choosing the correct 
equation in every practical case. We pay particular at- 
tention to the limits of the validity of known solutions 
and we carry out comparative measurements of magne- 
tostriction of the same thin film sample using both direct 
and inverse techniques. 



II. REVISION OF THE THEORY 

Already the existence of two "alternative" theories of 
the magnetostrictive effect creates some confusions. The 
anisotropic magnetostriction can be introduced in two 
ways. First, as the proper material deformations seen 
experimentally when the sample reaches the magnetic 
saturation (traditional technical definition 3 ). In the sim- 
plest case of an isotropic material are related only to 
the magnetization direction ( x axis in Fig. [1^,) and are 
given by a single magnetostrictive coefficient X s : 



+X S 



(1) 



In the second approach the magnetostriction is defined 
as the proper stresses originating from internal magnetic 
interactions (first principle approach!, see Fig. []})). For 
an isotropic material: 
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These definitions are equivalent providing that the in- 
ternal magnetostrictive stresses cr A - produce elastic defor- 



mations equal to the magnetostrictive ones eh, i.e. 



= A s 



E 



1 



(3) 



where v is the Poisson's ratio and E is the Young mod- 
ulus. 

In general, the initial magnetic state of the sample is 
not known and the only way to obtain the magnetostric- 
tion value is to saturate magnetization successively in 
two orthogonal directions (say x and y) and to measure 
the difference of deformation or stress between these two 
states (i.e. e x xx - e* y = 3/2 A s or a£ x - o* y = 3/ 2 cr A ). 

The corresponding energy density of interaction be- 
tween magnetostrictive deformation and elastic deforma- 
tion (magneto-elastic energy) can be obtained from two 
equivalent general expressions : w ms = — JT 



* 4 ° 



Very often the last expression is presented in a form 
where instead of a\ terms of type B 1 ' 2 cos 2 fa are used, 
where fa is the angle between the magnetization and the 
i-axis&i£. 

Combining all together for our isotropic case we obtain 
the angular variation of the magnetoelastic energy den- 
sity for magnetization rotation in the basal plane [x,y): 



Wrasifa) = - ^ K{a xx ~ <J yy ) COS 2 ((f>) 

3 

= (e xx - e yy )cos 2 {<t)) 

= B ia {e xx - Cyy) cos 2 (0) 



(4) 



where <f> is the angle between the magnetization and x 
axis. 

From this equation one can see that the physical mean- 
ing of B 1 ' 2 referred usually as the magnetoelastic cou- 
pling constant represents the characteristic magnetostric- 
tive stress <7 A illustrated in figure [T] B 1 ' 2 = — 3 /2 a x — 

For a thin magnetostrictive film deposited onto a 
substrate the planar components (e xx , e vy and e xy ) of 
the elastic deformation are fixed by the corresponding 
components of elastic deformations on the substrate 
surface. For a cantilever substrate elastically bent along 
its length x these elastic deformations are proportional 
to the curvature 1/R. A general expression of the energy 
of the unit area of such a sample in terms of 1 / R can be 
given as follows: 
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FIG. 1: Comparison of two definitions of the anisotropic 
magnetostriction, a) Spontaneous magnetostrictive deforma- 
tions related to the magnetization vector oriented along x 
axis. The relative deformation e^ x denned as the magne- 
tostriction constant \ s is accompanied by half smaller oppo- 
site transverse deformations Cy y and e zz . The case of positive 
A s when the initial cube of isotropic material is transformed 
to the elongated parallelepiped is shown, b) Alternatively the 
origin of magnetostriction can be defined trough the appear- 
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ance of internal stresses a xx — ^ , u yy 
which can be equilibrated by external stresses (big white ar- 
rows) (Tij = — cry to block the deformation of the material 
and, thus, to provide the elastic deformations opposite to the 
magnetostrictive ones. 



thickness) and the magnetic anisotropy energy density of 
the film is W an /t. 

For simplicity we consider only the uniaxial magnetic 
anisotropy with easy axis along the sample length. Here 
K is the usual anisotropy constant. 

By analogy with the bulk measurements the direct way 
to get the films magnetostriction is to turn the magne- 
tization direction from <f> — to 90° and to observe the 
corresponding change of the equilibrium substrate curva- 
ture that can be calculated by minimization of eq. ([5]) : 



R Ad 



(6) 
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cos 2 (0) - cos 2 (cV)(5) 



where W e i is the integral of the elastic energy density 
over the substrate thickness d; W ms — t w rns (t - thin film 



Practically this effect can be measured by deviation of 
a laser beam reflected from a point near to the free sample 
end: A8 = 2LA (-^). Here L is the sample length (more 
precisely the distance between the support and the laser 
spot on the sample). 

Alternative way to measure the magnetostriction is the 
more frequently used inverse technique: (-4) is modified 
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and according to eq. ([5]) the corresponding changes of 
the effective uniaxial anisotropy constant is measured: 



K eff =K + B± 



(7) 



As follows from eq. © and ([7]) both the direct and 
the inverse techniques should give the same result: the 
value of coefficient B where the magnetostrictive param- 
eters of the film are "hidden". The equivalence of these 
techniques is rarely discussed because they are practically 
never used together. 

A closer look to the problem shows that the parame- 
ters A and B may be different for different experiments 
and thus comparison of data obtained from different tech- 
niques is not so obvious. Moreover not only direct and 
inverse experiments should be compared with care but 
also data obtained from the same technique should be 
treated differently depending on the geometry of the sub- 
strate bending in the given experiment. 

In general, the Young moduli and Poisson's ratios of 
the substrate and of the film material are different and 
should be distinguished when used together. Below we 
demonstrate how different algebraic forms of the Pois- 
son's ratio of the cantilever substrate v s enters to A and 
B in two characteristic cases: anticlastic bending and 
cylindrical bending. And we show in a simple way why 
the Poisson's ratio of film Vj and Young moduli of both 
the substrate E s and the film Ef appear always in the 
same form for all kind of deformations. One can easily 
understand the variety of these forms of v {v s or vf) by 
reviewing characteristic cases of elementary deformations 
presented in Appendix [Al 

Small bending of a cantilever in one direction in- 
evitably induces an opposite bending in perpendicular 
direction in order to satisfy the absence of transversal 
stresses at its free side edges* 19 . In this case the principal 
surface deformation is followed by the transversal defor- 
mation as for an elastic element under uniaxial stress 
(see configuration 1 in Appendix [A"| . These two opposite 
bending produces an anticlastic shape of the cantilever 
that can not develop simultaneously with increasing de- 
formation. Due to the corresponding geometrical con- 
strains the process becomes non linear and a much more 
complicated deformation appears for large cantilever cur- 
vature. This effect is usually not considered in the papers 
on thin film magnetostriction while it can be non negli- 
gible as shown below. 

The conditions when bending can be considered small 
and the corresponding equations are linear can be found 
only in specialized papers on the mechanics of thin plates. 
The relevant critical parameter for thin cantilever with 
large length L to width W ratio is (3 ~ W/ V Rd. 

The linear transversal anticlastic deformation across 
the sample width z = y 2 u s /R is valid when this value 
at the sample edge y = W/2 is less than about QAd for 
v w 1/3 i.e. [3 < 1 (see Appendix IB"]) . 

At larger (3 the anticlastic profile cannot progress any- 
more in the same way because an additional edge axial 



tension created by it can no longer be neglected. The 
shape of the plate flattens and for (3 > 3 one can con- 
sider the bending as purely cylindrical. Another case of 
the cylindrical bending is when transverse deformation 
is blocked near to the sample fixed end by the support 
(clamping effect^). An immediate consequence of the 
fixed transversal deformations is the appearance of the 
transversal stresses a yy (configuration 2 in Appendix [A"]l . 

It is obvious that both the effective elasticity coeffi- 
cient of the plate bending A and the coupling coefficient 
B between the bending and the magnetostrictive stress 
in the film are different in these two limit cases. 

Pure bending in the principal cantilever direction x 
produces the longitudinal elastic deformations e xx (z) = 
z/R (see Fig. EJ. 

In the anticlastic case the elastic energy of the unit 
area of the curved plate is 



,1/2 



1 
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24 \R 



(8) 



taking into account the relation between e xx and a xx 
(configuration 1 in Appendix [Al a xx = Ee xx ). This gives 
A = 1/12 E s d 2 . 

The magnetoelastic coupling energy of the thin film 
with planar magnetization is fully determined by the sur- 
face deformation of the much thicker substrate: 



tB->' 2 I 



e w)z=d/2 C0S 



tW - 2 ^-(l + Vs ) cos 2 4> 
ZK 



(9) 



that gives B = 1/2 dB 1 - 2 (1 + v s ). Thus expression (JT]) 
of the effective magnetic anisotropy for inverse magne- 
tostriction measurements in the anticlastic case includes 
the Poisson's ratio of the substrate in the form (1 + v s ) . 

From eq. © we obtain the variation of the curvature 
l/R produced by the magnetization rotation from <f> — 
to 90° to be used in the direct technique: 



R 



anticlastic 



Bt 1 + v s t „ 9 , s 
Ad- 6 ^f^ B ^ ™ 



This result exactly corresponds to the solution given 
by Tremolet de Lacheisserie^. 

We should remember here that the Young's modulus 
E s and Poisson's ratio v s of the substrate enters into 
W e i and W ms while B 1 ' 2 includes the Young's modulus 
Ef and the Poisson's ratio vt of the film. 

For the cylindrical curvature, one have to use the re- 
lation between a xx and e xx corresponding to the fixed 
transversal deformation (e yy — 0, configuration 2 in Ap- 
pendix [X] cr xx = E/(l — v 2 ) e xx ). In this case, the in- 
tegration ([5J) gives harder plate bending coefficient A = 
1 /i2 Ed 2 /(l — v 2 ). At the same time the elastic coupling 
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becomes smaller B = 1/2 dB 1 ' 2 . Consequently, v s does 
not enter into the equation ([7]) of K e f t used for inverse 
measurements unlike the previous anticlastic case. 
As to the direct measurements, eq. © provides: 

This equation giving (1 — i/ s ) smaller substrate cur- 
vature than in the anticlastic case can be found in the 
literature, often without any justification. In reality, it 
has no practical application because in direct measure- 
ments the cantilever curvature is extremely small (/3 <C 1) 
and clamping effect can be neglected for cantilevers with 
large L/W ratio. 

One should note that the elastic deformation e xx due 
to the plate bending is not the only deformation produced 
by the magnetostrictive stress during direct measure- 
ments of magnetostriction. This stress also changes the 
sample length creating homogeneous deformation e xx . 
Consequently the position of the neutral line - where to- 
tal longitudinal deformation e xx {z) + e xx = - is shifted 
as discussed by many authors (see for example ^ 14 ' 20 . 
The elastic energy corresponding to this additional de- 
formation does not interfere with the energy of the pure 
substrate bending and, thus eq. © is unaffected. The 
shift of the neutral line can be illustrated by consider- 
ing separately symmetrical and antisymmetrical defor- 
mations produced by films deposited on both surfaces of 
substrates (see Appendix ICl) . This small e xx « X s t/d 
usually does not produce any measurable effect. 



netic field. The sample can be rotated in its plane in or- 
der to measure the film local magnetization either along 
its length (as shown in the figure) or along its width. 

The sample is mounted on a special bending jig shown 
Figl^l One end of the sample is glued on a support so 
that the film is on the top free surface from which the 
laser is reflected. The other sample end can be moved 
vertically (zl) by a screw in order to bend it up or down. 
To respect the sign given in the theoretical equations we 
define positive curvature when the sample is bent down- 
wards, and negative when the sample is bent upwards. 
We have found that the most precise way to get the sam- 
ple curvature in our experimental setup is to analyze the 
modification of the reflected laser beam. The light spot 
observed on a screen after reflection from curved surface 
is spread in two directions proportionally to the prin- 
cipal curvatures. This method is similar to the indus- 
trial multi-beam technique of measuring curved surfaces. 
With a single beam we can measure the curvature at 
the same small area where the MOKE hysteresis is stud- 
ied. In order to justify and calibrate our single beam 
technique we have analyzed the reflection from a number 
of cylindrical objects of known radius. Measuring our 
sample we could see that the light reflected near to the 
support forms on the screen a narrow line thus showing 
that the sample is curved here only along its length and 
remains flat in the transverse direction (cylindrical cur- 
vature). Further from the support we could see that the 
laser spot becomes elliptical thus demonstrating the an- 
ticlastic shape of the sample deformation (see Appendix 

m. 



III. EXPERIMENTAL DETAILS 

We apply both direct and inverse techniques to 
the same sample. For illustration we have chosen a 
{[Tb 34% Co 66% ]60 A/[Co 42% Fe 58 o /o ]50A}x 10 multilayer 
deposited onto rectangular Corning glass substrate (22 x 
5 x 0.145 mm 3 ) from CoFe and TbCo mosaic 4 inch tar- 
gets using a Z550 Leybold RF sputtering equipment with 
a rotary table technique. An uniaxial anisotropy along 
sample length was induced by the field of a permanent 
magnet applied during deposition process. The direct 
measurements were realized using cantilever deflectome- 
try technique. Details of the sample preparation and of 
the cantilever measurements can be found in 23 . 

For the inverse measurement we have studied local 
magnetization curves of bent sample using the magneto- 
optical Kerr effect (MOKE) as shown in figure O We 
have chosen the longitudinal Kerr effect geometry for 
which the planar magnetization component parallel to 
the plane of incidence is measured. The incident laser 
light is polarized vertically and the light reflected from 
the film surface is split by a Wollaston prism to 2 pho- 
todetectors so that their signal difference is proportional 
to the Kerr rotation. The incidence plane and the sample 
surface are parallel to the direction of the applied mag- 




Polarizer 
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H 



FIG. 2: (Color on line) Schematic diagram of the indirect 
measurements of the film magnetostriction using MOKE mag- 
netometer. The geometry of sample bending is shown at the 
top of the figure. The dashed line shows the neutral line where 
the longitudinal deformation is zero. 
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IV. RESULTS AND DISCUSSION 
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FIG. 3: (Color on line) MOKE magnetization loops mea- 
sured close to the sample support (x — 2 mm) for different 
positive (left column) and negative (right column) cantilever 
curvatures. Red full square symbols - magnetic field applied 
along the sample. Blue open circle symbols - magnetic field 
applied across the sample as shown in c) and c') near to the 
corresponding curves 



For positive curvatures (left column), the large axis 
hysteresis loop shape does not vary, keeping the rect- 
angular shape characteristic for the easy axis magneti- 
zation, whereas the short axis hysteresis loop is getting 
more and more slanted when curvature increases. Thus 
we clearly see that a positive curvature (stress) reinforce 
the anisotropy. Negative curvature produces the opposite 
effect: the uniaxial anisotropy, first, is reduced by the ap- 
plied stress (the slope of the short axis hysteresis loop is 
reduced by stress counterbalancing the initial anisotropy) 
and, then, the easy anisotropy axis switches from longitu- 
dinal to transversal direction and the anisotropy increases 
again. In our sample we observe that this switching is not 
abrupt: the long axis hysteresis starts to slant before the 
short axis hysteresis becomes completely rectangular (see 
Fig. [3] (b')). As we show below, this effect appears if the 
sample is not perfectly homogeneous. 
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FIG. 4: (Color on line) Dependence of anisotropy field mea- 
sured as a function of curvature for fixed position [x = 2 mm) . 
Positive values of Hk (open blue circles are obtained for 
H || x) and negative values of Hk (full red squares are ob- 
tained for H || y). 



We have studied the variation of the MOKE hystere- 
sis with curvature at x — 2 mm from the support where 
the maximum curvature is achieved. Characteristic re- 
sults from a sequence of measurements done for different 
bending in positive and negative directions are shown in 
Fig. [3l For every curvature the sample has been ro- 
tated so as to obtain magnetization cycles both for the 
magnetic field applied along the sample large axis (full 
red square symbol) and perpendicularly to it (open blue 
circle symbol). 

As observed with VSM similarly to our previous 
study 23 the unstressed sample has a uniaxial anisotropy: 
rectangular magnetization along the sample length shows 
the easy axis and slanted magnetization along the sam- 
ple width (hard axis magnetization) provides the initial 
anisotropy value (see Fig. OUa)) 



Contrary to what was observed on the VSM cycles, 
back and forth branches in the hard axis loop are not fully 
superimposed and the cycle is somewhat asymmetrical. 
This is not an unusual feature of MOKE measurements 
caused by a non linear interference with the transverse 
magnetization signal 2 ^. We have found that for magne- 
tostrictive films this MOKE loop feature may have an 
additional reason. We observed a superposition of the 
magnetization and magnetostrictive response when the 
free sample end did not touch the screw at all and, thus, 
could move freely under the magnetostrictive stress in 
the film. Probably a small part of the latter effect can 
persist even when the sample end is touching the screw. 

This small distortion of the hard axis magnetization 
loop can be excluded by the following procedure of 
the determination of the uniaxial anisotropy field Hk- 



6 



The loop is firstly "symmetrized" by addition of the 
back and forth branches. Then the central part of it 
[-20 Oe; + 20 Oe] is linearly fitted and H K is de- 
termined as the intersection of the fitted straight line 
with the saturation magnetization level. The anisotropy 
field Hk and the corresponding anisotropy constant K 
(Hk = 2K/{hqM s ) where M s is the saturation magne- 
tization) are defined positive when the hard axis loop is 
seen for applied field parallel to the small sample axis (see 
Fig. [3]a)-c), and negative if the hard axis loop is seen for 
applied field parallel to the large sample axis (see Fig. [3] 
c')- 

For a perfectly homogeneous sample with a well de- 
fined uniaxial intrinsic anisotropy direction (i.e. long 
sample axis), we would observe an abrupt reversal be- 
tween slanted and rectangular loops at the moment of 
90° rotation of the easy axis when the initial anisotropy 
is compensated by the stress induced anisotropy. In re- 
ality, as we indicated above (Fig. Wp) the transition is 
gradual. There is an interval around the compensation 
point where Hk can be formally measured for both direc- 
tions. This is due to small variations of both magnitude 
and direction of the intrinsic anisotropy in the sample 
area probed by laser. This kind of anisotropy inhomo- 
geneity was already observed and quantified in our previ- 
ous work2£. There we have shown that large domains of 
opposite magnetization exist where maximum easy axis 
distribution is about ±6° from the average direction over 
the whole sample. In the intermediate region every mea- 
sured Hk corresponds not only to the average anisotropy 
of given sign but also to its relative occurrence. 

Figure [4] shows the evolution of the obtained positive 
and negative Hk with the substrate curvature (R^ 1 ) or 
the equivalent surface strain e xx — d/2R). 

In order to get the expected linear variation of Hk with 
curvature one has to consider only points outside the in- 
termediate compensation region shaded gray in the fig- 
ure. Inside this region the slope is about twice smaller 
because the repartition between the areas with positive 
and negative Hk also varies with curvature. So if this re- 
gion is not identified the resulting magnetostriction value 
can be considerably reduced. The slope of the line fitting 
the points outside the compensation region =8.2 
Oe.m = 650 A allows us to determine the magnetoelastic 
coupling coefficient of the film i? 7,2 . 

As mentioned above for this measurement the defor- 
mation is cylindrical close to the support. So, for eq. ([7]) 
we have to take the corresponding form of the effective 
magnetoelastic coefficient B: 

2 ij, M s dH K , v 

B =—wm ( } 

With HqM s — 1.4 T obtained from the magnetic mea- 
surements and d = 145 urn we found B 7 ' 2 = 6.3 MPa 

Far from the support the deformation of the sample 
gets an anticlastic shape if the bending is not too strong 



(see Appendix IB")) . As was mentioned above, we observe 
this anticlastic deformation as the dispersion of the re- 
flected laser beam: the laser spot is spread to an ellipse 
after reflection from the sample. For the deformation be- 
ing equivalent to the deformation of narrow cantilevers 
the magnetoelastic coefficient B is by (1 + v s ) larger. 
Correspondingly one have to use for this case another 
expression relating B 1,2 and the induced anisotropy: 

2 _ _}_IIoM 1 _3Hk_ 

~ 1 + v s d 8(1/ R) [ ' 

In order to observe the effect of this shape modifica- 
tion we have carried out measurements at different points 
from the support. 

MOKE hysteresis loops obtained for a given deforma- 
tion zl = +1.6 mm, left column, and zl = —1.6 mm 
right column, when laser spot is swept along the sam- 
ple length from a position close to the clamped tip 
(x = 2 mm) to x — 14 mm trough x = 8 mm are pre- 
sented Fig. [5j 

For positive deviation [zl = +1.6 mm), the closer to 
the clamped sample tip [x = 2 mm) the more anisotropic 
the sample is: The small axis hysteresis loop is more 
slanted when recorded closer to the clamped tip. The 
hard axis (small side of the sample) is thus getting harder 
when x decreases. Hk as defined above remains pos- 
itive decreasing with x from its maximal value to its 
initial value close to the free end. The reason for this 
is that the created curvature and the corresponding in- 
duced anisotropy reduce progressively with the distance 
x from the support completely disappearing at the free 
end x = L (see Appendix [C| . 

For negative deviation (zl = —1.6 mm) we observe 
the inversion of easy and hard direction at small x where 
the induced anisotropy overcomes the initial one. Cor- 
respondingly, Hk varies from a negative to a positive 
value along the sample length. This transition is not 
abrupt similarly to the modification of Hk as function 
of the bending curvature discussed above and there is an 
intermediate compensation area approximatively in the 
middle of the sample. We represent the evolution of the 
induced anisotropy along the sample length by the differ- 
ence between Hk obtained for positive (H^) and nega- 
tive {H K ) bending (see Fig. [B]). Such presentation allows 
us to exclude the inhomogeneity of the initial anisotropy 
of the undeformed sample. We can distinguish two kinds 
of points: when H K and H K have different signs in the 
part of the sample close to the support (open circles in 
the figure) and when H K and H K remain both positive 
in the second part of the sample (closed square in the 
figure). Similarly to the evolution with l/R discussed 
before these two kinds of points coexist in the interme- 
diate compensation area and expected theoretical value 
should lie somewhere between them. Except of the in- 
convenience caused by the sample inhomogeneity, there 
are two effects to be considered for quantitative analysis 
of these data: the effect of the support clamping and the 
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edge blocking for large curvature. 

The cylindrical shape imposed by the clamping at the 
support extends to some distance thus producing in this 
region the strengthening of its bending and the reduction 
of the induced anisotropy as discussed above (see discus- 
sion before eq. (fTTj) . The total shape of the clamped 
sample does not depend on the bending amplitude until 
the curvature is small enough and the linear equations 
can be used. The variation of the induced anisotropy 
due to clamping effect can be illustrated by eq. (|C9p of 
Appendix 
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this linear approximation (3 ~ W/V Rd < 1 is valid, 
for our sample is (1/R) m ax = 6 m _1 . In our exper- 
iment the largest sample curvature considerably over- 
comes this value near to the support (see Fig. [3]) and 
a further gradual modification of the anisotropy with 
curvature should be added according to the calculated 
correction k(1/R) (see Appendix [B] eq. (|B7j) and Fig. 
|Ht>))- Since 1/R reduces almost linearly along the sam- 
ple length, the low curvature approximation can be used 
till maximal allowed bending only close to the sample 
free end (x > 14 mm) where, unfortunately, the induced 
anisotropy is too small for precise measurements of the 
magnetostriction. 

Thus, in the main part of the sample we have a com- 
plicated superposition of at least three effects: the com- 
pensation effect with the revealed inhomogeneous initial 
anisotropy, the clamping effect and the non linear shape 
effect. 

Interestingly, our experimental measurements near to 
the cantilever support discussed above (Fig. [3J H]) have 
no such complication. The clamping of the cantilever end 
imposes the cylindrical curvature at any width to length 
ratio and at any bending amplitude, and excludes all 
other effects considered above. Consequently the deter- 
mination of the magnetoelastic coefficient £? 7 ' 2 has less 
uncertainty. 
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FIG. 5: (Color on line) MOKE magnetization loops measured 
for opposite sample bending zl = +1.6 mm (left column) 
and zl = —1.6 mm (right column) at different position on 
the sample: x = 2 mm, x = 8 mm and x — 14 mm. Red 
full square symbols - magnetic field applied along the sample 
length and blue open circle symbols - magnetic field applied 
across the sample 



FIG. 6: (Color on line) Modification of the induced anisotropy 
along the sample length. Open circles shows results involv- 
ing anisotropy field > created by positive bending and 
Hft < created by negative bending. Results shown by full 
squares involve positive anisotropy fields and for both 
opposite bending. A theoretical estimation taking into ac- 
count the clamping effect is shown by the dashed line. 



An example of such an estimation is shown by the 
dashed curve in Fig. [6] calculated with the following pa- 
rameters of our sample: sample length L — 18.5 mm, 
v s = 0.2, characteristic length of the clamping ef- 
fect a as W/2 = 2.5 mm and amplitude coefficient 
6FB~*' 2 /E s Wd 2 = 1.2 xl0 6 A/m 2 . The theoretical pre- 
dicted deviation from the linear dependence is clearly 
visible close to the sample support. 

The maximum curvature, for which the condition of 



For our direct measurements of magnetostriction by 
deflectometry the anticlastic shape of the sample can 
be always considered: The aspect ratio L/W ~ 4 of 
the cantilever is sufficiently large to neglect the effect of 
clamping^ and the curvature is extremely small. Thus 
corresponding anticlastic parameter j3 < 0.01 and only 
eq. (fTTj)) have to be used. 

Figure [7] shows deflectometry measurements of our 
sample under the magnetic field applied in two directions: 
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parallel to the sample length (easy axis) and perpen- 
dicularly to it. The easy axis magnetization, as should 
be, does not produce any magnetostrictive effect. The 
maximum laser deflection in the saturating perpendic- 
ular field is AO — 0.15 mrad giving the sample curva- 
ture A(l/R) = A6/2L = 4.3 xfCTV- 1 (distance be- 
tween the support and the laser spot on the sample is 
L = 17.5 mm). From these measurements we obtain 
B 1 ' 2 = 6.8 MPa using eq. (| 10[) with following param- 
eters: E s = 60 GPa, v s = 0.2, d = 145 |im, t = 110 nm 
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FIG. 7: (Color on line) Deflectometry measurements of 
the magnetostrictive bending under the field applied along 
the sample length and perpendicularly to it as schematically 
shown near each corresponding curve. The geometry of the 
laser beam deflection is sketched in the inset. 

This value is in good agreement with the value ob- 
tained above (6.3 MPa ) using the inverse technique, ft 
is difficult to get better correspondence of two techniques 
taking into account that there are parameters which enter 
to the calculation for only one of the two techniques: for 
the direct technique - elastic moduli of the substrate E s 
and v s , sample length L, other dimensions of the instal- 
lation and the film thickness f; for the inverse technique 
- saturation magnetization of the film M s , curvature cal- 
ibration 1/R and determination of Hk- Another impor- 
tant difference is that the substrate thickness enters as 
d 2 to the direct technique and as 1/d to the inverse one. 
Eventual influence of the inhomogeneity on B 1 ' 2 also can 
not be excluded: the direct technique provides the av- 
erage value over the whole sample whereas the inverse 
technique measures a local value. 



experimental conditions are correctly taken into account. 

It is natural to characterize the thin films magne- 
tostriction by <7 A since the lateral deformations of the 
film are blocked by the substrate and one can measure 
only the stress (see Fig. [TJd). This is in contrast to mea- 
surements of the magnetostrictive deformation \ s of bulk 
materials (see Fig. [T]j). 

The key factor to determine the exact expression for a x 
to be used is the shape of the curved cantilever realized 
in the particular conditions of the experiment. Both the 
elastic and the magnetoelastic energies depend on the 
peculiarities of the curved cantilever shape. Normally a 
small cantilever bending involves anticlastic motion of its 
side edges. 

This motion can be hindered by two effects: clamping 
by the support and non linear geometrical edge block- 
ing for large cantilever bending. The critical curvature 
above which this cylindrical bending with completely flat 
transversal profile appears is W 2 /Rd > 10 (more pre- 
cisely anticlastic parameter [3 > 3). For the direct mea- 
surements the critical curvature can be achieved only for 
extremely thin cantilevers with d < ^/X s W 2 t. This crit- 
ical value is obtained from eq. (|10|) neglecting the dif- 
ference of the elastic moduli of the film and the sub- 
strate. For A s = 10~ 4 , W = 5 mm and t — 110 nm, as 
an example, this gives d < 6 u\m. The indirect measure- 
ments are limited by the maximum strain that can be 
obtained on the cantilever surface without risk to break 
it (e m ax ~ 2 x 10 -3 ). So, in this case the critical limit to 
the cantilever thickness is less severe: d < \/0.2e m axW 
i.e. d < 100 |J.m for the same numerical values. 

In the extreme case of strongly curved thin cantilever 
clamped at one end the result of the local measurements 
of the induced anisotropy will vary not only with the dis- 
tance from the support but also across the sample width: 
K = K{x,y) (Appendix O 

Comparing the advantages of the two techniques one 
can think that the indirect technique is more sensitive 
and easier to realize and practically unlimited by the film 
thickness. This is true only for the materials with weak 
magnetic crystalline anisotropy. When this anisotropy 
is much larger than the maximum induced anisotropy 
\ s Efe ma x this technique loses its advantages. The mag- 
netostriction measurement of highly anisotropic materi- 
als by the direct technique is limited mainly by the avail- 
ability of strong enough magnetic field (H > Hk). In or- 
der to measure materials with low X s or very thin films by 
the direct technique one can reduce the substrate thick- 
ness. 



V. CONCLUSION 

The comparison of direct and inverse (" indirect" ) tech- 
niques of the film magnetostriction measurements shows 
that both techniques give exactly the same magnetostric- 
tive stress <7 A (or B 1 - 2 = — 2/3<t a ) provided all different 



Acknowledgments 

D. Dekadjevi is acknowledged for the development of 
the MOKE experiment. 



9 



APPENDICES 

APPENDIX A: EFFECT'S OF THE POISSON'S 
RATIO 



potifi cmi'flfinii 1 


Free edge load: 

definition of the 

Young's modulus E 

and the Poisson ratio v 
p _ 

c xx — £ 

^yy — v&xx 




Oyy = 


configuration 2 


Fixed edge load 
e - 

V XX 

^yy — ~\~v(Jxx 




p - 


configuration 3 

/ A ^xx 


Antisymmetric 
load / deformation 

P — —p — l+lJ n — 

CXX Cyy E OxX 2 Q 

(G- shear modulus) 




configuration 4 


Symmetric load/deformation 

&xx — ^yy — ~£j~o~xx 

(Case of isotropic lateral 

stresses) 




U Oyy = + 0^ 



TABLE I: Summary of characteristic appearances of Poisson's 
ratio v 



Even in the simple case of isotropic materials the ap- 
pearance of the Poisson's ratio v in different formula re- 
lated to the magneto-elastic interactions presents some- 
times a " mystery" . In order to simplify this issue we give 
here a simple summary of relations between stresses cr^ 
and deformations = x, y, z) of an isotropic mate- 

rial for typical conditions of mechanical loads useful for 
the thin film geometry with their proper characteristic 
appearance of v. Let's consider a small square element 
with the principal load/deformation along the horizontal 
x axis and z axis perpendicular to the film surface where 
the normal stress u zz = 0. 

When the side edges of the element are free (a yy = 0) 
the lateral linear deformations e xx and e yy under uni- 
axial stress a xx define the Young's modulus E and the 
Poisson's ratio v (see configuration 1 in table [A"|l . The 
ratio e yy /e xx — — v is directly related to the anticlastic 
bending of narrow cantilevers with free side edges. 

The elastic properties of the element with fixed side 
ed ges e yy — can be obtained through a linear super- 
position of two crossed loads (o~ xx and <J yy ) with free 
edge conditions presented above, where o yy compensates 



the transverse deformation produced by principal load 
o~ xx '■ o- yy = vo~ xx (configuration 2 in table . The re- 
sulting strengthening u xx /e xx of the axial deformation 
(1 — v 2 )^ 1 corresponds to the strengthening of strongly 
bent wide cantilevers where the condition e yy — is due 
to the transverse flattening caused by the geometrical 
edge blocking (see Appendix [B| . 

Another useful superposition of two uniaxial loads a xx 
and a yy — ~(J XX (antisymmetric case) is equivalent to the 
pure shear deformation along the diagonal of the clement 
(configuration 3 in table [X]) . It appears in the magne- 
toelastic coupling where only the difference between two 
states with stress axes rotating by 90° is essential. 

Consequently the elastic parameters of the material 
appear in the form E/(l + v) in these cases (see eq.© 
and (fT0f ). 

For completeness we present also the symmetric super- 
position of o~ xx and a yy — a xx which correspond to the 
thermal and epitaxial stresses of thin films (configuration 
4 in table . Factor (1 — v) obtained here instead of the 
factor (1 + v) in the previous antisymmetric configura- 
tion enter to the expression of the substrate curvature 
(eq. pop ) for this case well known as Stoney's formula 15 . 



APPENDIX B: EVOLUTION OF ANTICLASTIC 
DEFORMATION 

The details of induced anisotropy of magnetostrictive 
films on curved substrates related to their anticlastic de- 
formation, clamping effect were largely discussed in the 
literatur o 6 ' 9 ' 14 i 20 . All these numerical and analytical so- 
lutions were limited to very small deformations. It is 
interesting to consider the effect of the evolution of the 
form of a cantilever with increasing curvature. The an- 
ticlastic deformation can not be maintained for strong 
curvatures because of appearance of tear stresses at the 
side edges of the cantilever. General equations describ- 
ing thin plate 2D profile Z(x, y) are non-linear. For our 
cantilever with principal curvature 1/R along its length 
(Z(x, y) = — x 2 /2R+z(yj) there is an analytical solution: 
the transverse profile of the vertical deviation z{y) can 
be obtained from the 4th order differential equation^: 

where 

/3 2 = ^[3(1-- 2 )] 1/2 (B2) 

is a parameter (anticlastic parameter) which determines 
the degree of the non-linearity . 
W is the cantilever width and d its thickness. 
The solution of this equation is : 
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width and its amplitude saturates (Fig. [9]-a). This edge 
deformation region is clearly demonstrated on images of 
bent metallic sheets shown by Conway and Nickola 4 . The 
transition from anticlastic to cylindrical bending pro- 
duces corresponding modification of the ratio between 
induced anisotropy and the principal cantilever curvature 
l/R [See changes of the forms of the effective coefficient 
of magnetostrictive coupling B in eq. (fit)]) and (TTTj) ]. 

It is interesting to follow the whole evolution of this ra- 
tio from small (linear regime) to large (non linear regime) 
cantilever bending. 

The induced anisotropy is determined by (e xx — e yy ) 
at the cantilever surface where the film is deposited: 



!JV 



dd 2 Z _ dl 
2lh? ~ ~2R 
dd 2 Z _ dd 2 z 
2V ~ 2d^2 



(B5) 
(B6) 



For a narrow cantilever e yy = —ve xx is valid for any 
curvature ((3 <C 1 see eq. ([B4I0 . 

The induced anisotropy measured in the center of wide 
cantilever of the same curvature must be corrected by a 
factor: 



FIG. 8: (Color on line) Anticlastic deformation of a bent can- 
tilever for different values of the reduced curvature W 2 /Rd. 
a) small enough bending W 2 /Rd — 1 (anticlastic parameter 
13 = 1.28) b) strong bending W 2 /Rd = 10 (/3 = 4) with a flat- 
tening in the center. Transverse profile calculated from eq. 
(|B3|) . (z/d) is differently magnified in a) and b) for clearness. 
We added to this figure a gradual flattening qualitatively re- 
producing the clamping effect of the support (at x = 0). 



^^{G- cosh cos (/?#) 
+C + sinh(/?^)sin (/?#)} 



(B3) 



with 



C- 



sinh(/3/2) cos(/?/2) =p cosh(/3/2) sin(/?/2) 
sinh(/3) + sin(/3) 



For j3 <C 1, it gives the anticlastic transversal curva- 
ture (v I ' R) opposite to the principal curvature of the can- 
tilever (-1/-R): 



W 2 /8) 



(B4) 



as illustrated in Fig. [8] It is interesting to note that this 
parabolic transversal deformation proportional to l/R 
persists till (3 ~ 1. 

For larger (3 the cantilever starts to flatten in the cen- 
ter. For (3 > 3 (see Fig. [5] ) the anticlastic deformation 
remains only at the sample edges in a region of W/(3 



n(l/R) = 



(1 + v)e xx 



1 



„ d 2 z 
1 + R ^ 



(B7) 



V=0j 



For v = 1/3 the evolution of k is shown in Fig. [5]- b). 

When bending is small (l/R < W/d 2 or < 1 ) the 
anticlastic shape does not depend on the cantilever width 
d 2 z/dy 2 = v/R and k = 1. When bending is large 
(l/R > 10W/d 2 or (3 > 3) and the center of the cantilever 
becomes flat (d 2 z/dy 2 « 0) we have to correct the usual 
expression of the induced anisotropy by k = 1/(1 + v). 



APPENDIX C: STRESS VARIATION BY THE 
CANTILEVER THICKNESS AND LENGTH 

There is a large variety of the stress distribution for 
different cases of the cantilever deformation. We con- 
sider here two specific cases relevant to our experiment. 
In this appendix, for simplicity, we discuss only the free 
edge linear curvature of the cantilever corresponding to 
configuration 1 in Appendix |A1 with a xx = E s e xx . Here 
we have to distinguish Young modulus and Poisson's ra- 
tio for the substrate (E s and v s ) and for the magne- 
tostrictive film (Ef and Vf). The latter appear only in 
B 1 ' 2 = -3/2 A s E f /(1 + VJ ). 

The deformation of the cantilever under the magne- 
tostrictive stress of the film with magnetostriction +A S 
deposited onto only one surface is not symmetrical with 
respect to its midplane. The corresponding effect of the 
shift of the neutral line of the elastic plate deformation 
somewhat complicates the analysis provided in different 
papers*^. This practical problem can be considered as 
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FIG. 9: (Color on line) Effects of the reduced curvature w2 /nd: 
a) Transverse deviation z at the edge of the cantilever (y — 
W/2) relatively to its thickness d. b) Correction for induced 
anisotropy in the center of a wide cantilever compared to an 
equally bended narrow cantilever, v = 1/3 is used for the 
calculations. 



the linear superposition of 2 simpler problems: symmet- 
ric - the deformation of the substrate by 2 equivalent films 
deposited on both surfaces with magnetostriction coeffi- 
cient +1/2 A s , and antisymmetric - 2 films with opposite 
magnetostrictions coefficients +1/2 A s and —1/2 A s on 
opposite cantilever surfaces (see Fig. [TO] ) . 

Let's consider the influence of the longitudinal stress 
aix = v X of the film magnetized along the x axis as 
indicated in Fig. [TO] 

In the symmetric case, evidently, the only deformation 
is lengthening e xx and there is no curvature. From the 
balance of the forces between the constant magnetostric- 
tive stress on both surfaces and the Hooke's force of the 
sample lengthening one gets: 



t <j' 



dE„ 



(CI) 



In the antisymmetric case the pure curvature (1/i?) 
with the neutral line at the middle of cantilever thickness 
is produced without any lengthening of the sample, i.e. 
e-xx(z) = z/R. Here one should balance the momentum 
of surface magnetostrictive forces and the momentum of 
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FIG. 10: The practical experimental configuration with a 
film deposited on one substrate surface can be considered as 
the superposition of a symmetric problem with 2 equivalent 
films deposited on both surfaces (sample lengthening only) 
and an antisymmetric problem with 2 films of opposite mag- 
netostriction on both surfaces (sample curvature only). 



the bulk Hooke's forces relative to the neutral central 
line: 



fd/2 

t a x d/2 = / z E, 

J-d/2 



c (z)dz 



that gives the resulting curvature: 
1 Qta x 



R(aL) 



(C2) 



(C3) 



For the initial combined problem, the deformation of 
curvature and lengthening have the same origin and are 
related by the following relation: 



6 (£) £ 



(C4) 



Thus the maximal value of e xx (z) is 3 1 e xx \ and the neutral 
line of the total deformation e xx (z) + e xx is located at 
z = -d/6 (See Fig. [TO]). 

For complete description of the effect of the magne- 
tization rotation more components of the stress should 
be considered: the tranversal stress for the magnetiza- 
tion along x axis fr| y = — 1/2 a x and the stresses for 
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the perpendicular magnetization <r^ x 



1 /2 cr A and 



cr . The position of the neutral plane for each 
component is the same and their total influence increases 
curvature l/R(al x ) by 3(1 + v s )j2 found in eq. (fTU)) . 
The shift of the neutral plane is found not only for mag- 
netostrictive film stres o 9 ' 14 i 20 i 24 but also for thermal or 
epitaxial stresses (— and references therein). 

Second question we have to address is the stress distri- 
bution and the corresponding curvature variation along 
the cantilever length. The simplest way to obtain this is 
to balance the momentum of forces around a given point 
of the cantilever curved by force F applied to its free end 
(x = L), the other end (x = 0) being clamped. 

Let's consider the equilibrium of part [x, L] of the can- 
tilever around the middle point at its left edge. The 
momentum of bulk Hooke's forces (see the second term 
of eq. (|C2|) ) applied to x should be equal to the applied 
force momentum. 



d/2 p 

z E s e xx (z)dz = —(L — x) 

-d/2 



w 



that gives: 



1 



12F 



R E s Wd 3 



(L-x) 



(C5) 



(C6) 



This simple solution is valid when the sample side 
edges are free and the anticlastic deformation occurs ev- 
erywhere along the sample length. In reality the anticlas- 
tic deformation is blocked at the sample support (clamp- 
ing) with corresponding bending hardening by a factor 



1 — v 2 ( see configuration 2 in Appendix [Al and modifica- 
tion of elastic coefficient A in eq. (O from anticlastic to 
cylindrical bending). This phenomena can be described 
qualitatively by introducing an effective Poisson's ratio 
correction v{x) = v s exp(—x/a) where a ~ W/2 is the 
characteristic relaxation length of the clamping effect. 
With this correction: 



1 



12F 



R E s Wd 3 



(l-i> 2 (x))(L-x) 



(C7) 



The clamping produces even stronger influence onto 
the induced magnetic anisotropy that has to be multi- 
plied by (1 + v) according to modification of the effective 
magnetoelastic coefficient B in eq. (O from anticlastic 
to cylindrical bending. So we obtain: 



V ; 2 l + v{x) 



(C8) 



(see eq. ©). 

Combining 1/R(x) and B{x) in eq. ([7]) we get a non- 
linear variation of the anisotropy with the distance from 
the cantilever support: 



QFB~<> 2 



K eff (x)=K + 
{l + v s ){l-v{x)){L-x) 



E P w*-:r — 7 . (C9) 

instead of the traditionally considered linear variation 
without clamping effect [y = 0). 
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